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In this work we perform a comparison between Classical Molecular Static (CMS) and quantum 
Density Functional Theory (DFT) calculations in order to obtain the diffusion coefhcients for diluted 
Fe-Cr alloys. We show that, in accordance with Bohr’s correspondence principle, as the size of 
the atomic cell (total number of atoms) is increased, quantum results with DFT approach to the 
classical ones obtained with CMS. Quantum coherence effects play a crucial role in the difference 
arising between CMS and DFT calculations. Also, thermal contact with the environment destroys 
quantum coherent effects making the classical behavior to emerge. Indeed, CMS calculations are in 
good agreement with available experimental data. We claim that, the atomic diffusion process in 
metals is a classical phenomena. Then, if reliable semi empirical potentials are available, a classical 
treatment of the atomic transport in metals is much convenient than DFT. 
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To characterize the crossover between the quantum 
and the classical worlds is a fundamental quest of 
modern physics. For some systems, it is clear that, 
classical physics arises from quantum physics in the 
large-number limit. This is the Bohr’s correspondence 
principle [l|. 

It is not fully understood yet how the many-particle 
limit gives rise to classical physics, and how much of 
quantum physics still remains. In this framework, the 
environment and its temperature plays an important 
role for the decoherence process (2|. However, quan¬ 
tum coherence effects at ambient temperatures where 
observed in biological systems Q. 

Also, for many body systems we must deal with 
both electrons and atomic nuclei dynamics. Due to 
their masses the nuclei move much slower than the 
electrons. Then in the Born-Oppenheimer approxi¬ 
mation the nuclei generate a static external poten¬ 
tial in which the electrons are moving. 

First principles {Ab initio) quantum mechanical 
method, such as DFT, are employed to obtain the 
electronic behavior. DFT reduces the quantum many- 
body problem to the use of functionals of the electron 
density [^ , it is presently the most successful approach 
to compute the electronic structure of matter. Its ap¬ 
plicability ranges from atoms, molecules and solids to 
nuclei and liquids @ . The electron density determines 
the potential energy surface that represents the force 
field where the nuclei dynamics occurs. 

In a different approach, a CMS treatment, employs 
phenomenological semi-empirical potentials in order 
to estimate this force field. Then, the quantum na¬ 
ture of electronic structure is not taken into account. 
There is a wide variety of semi-empirical potentials, 
which vary according to the atoms being modeled. Ab 
Initio simulations take into account the quantum na¬ 
ture of the electrons, which implies in a higher compu¬ 
tational cost than CMS. Hence ab initio simulations 
are limited to smaller systems. 

In this work atomic diffusion, in Fe-Cr diluted al¬ 
loys is studied with both CMS and quantum DFT cal¬ 
culations, in the context of a multi-frequency model. 
We show that, in accordance with Bohr’s correspon¬ 
dence principle, as the total number of atoms is in¬ 
creased, quantum results with DFT recover the clas¬ 


sical ones obtained with CMS. 

Quantum coherence effects, which are only taken 
into account by DFT, play a crucial role for both the 
convergence issue of the DFT results as a function of 
supercell size, as well as, in the difference arising be¬ 
tween classical CMS and quantum DFT calculations. 

In addition, CMS calculations are in good agree¬ 
ment with available experimental data for both solute 
and solvent diffusion coefficients. This may not be 
surprising for a macroscopic system especially for high 
temperatures that destroys any quantum coherent ef¬ 
fect. 

Diffusion plays an important role in the kinetics 
of many materials processes. Experimental measure¬ 
ments of diffusion coefficients are expensive, difficult 
and in some cases nearly impossible. A complimen¬ 
tary approach is to determine diffusivities in materi¬ 
als by atomistic computer simulations. In addition to 
predicting diffusion coefficients, computer simulations 
can provide insights into atomic mechanisms of diffu¬ 
sion processes, creating a fundamental framework for 
materials design strategies. 

Also, Fe-Cr alloys at low temperatures has impor¬ 
tant technological consequences. Due to their good 
resistance to void swelling 11 , 0 , Fe-Cr based alloys 
are of special interest for nuclear applications (in Gen¬ 
eration IV and fusion reactors). 

Atomic transport theory allows to express the dif¬ 
fusion coefficients in terms of the atomic frequency 
jumps, this is commonly known as the multi-frequency 
model 0 . Recently, attends were made in order to de¬ 
scribe the diffusion process by obtaining numerically 
the needed frequency jumps with DFT calculations. 
Although disagreement between the experimental and 
ad-initio based calculated diffusion coefficients where 
observed in bcc alloys such as Ni-Cr and Ni-Fe 
and for aFe-Ni and aFe-Cr alloys K| as well as for 
Mg, Si and Cu diluted in fee Al ll|. However, in 
a recent work Huang et al. [I^, for Fe based di¬ 
luted alloys, have performed DFT based calculations 
for the tracer diffusion coefficient with a larger num¬ 
ber of atoms (128 instead of 54 as in 0) that are in 
good agreement with the experimental data. 

On the other hand, one of us has recently shown 
[l^ . that tracer diffusion coefficients performed with 
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CMS based calculations in diluted Ni-Al and Al-U 
fee alloys are in excellent agreement with available 
experimental data for both systems. 

We focus here on the tracer self- and solute diffu¬ 
sion coefficients in a binary ^4-5 alloy in the diluted 
limit. For diffusion mediated by vacancies analytical 
expressions, in terms of the frequency mmps, where 
calculated by Allnatt and Le Claire H for foe and 
bcc lattices respectively. In the 2nd-nearest-neighbor 
binding model, we identify the jumps as in Fig. [T] 
The self-diffusion coefficient can be written as, 

D\ = a^uJoCvfo, (1) 


where a is the lattice parameter, ojq is the atom- 
vacancy exchange frequency in pure A and /o, the 
self-diffusion correlation factor that is fo = 0.7272 Q 
or fo = 0.7814 0, for bcc or fee metals respectively. 
At thermodynamic equilibrium the vacancy concen¬ 
tration depends on the temperature T as, 

f eY - TSY\ 

O'' = ) <"> 


where ks is the Boltzmann constant while Ej and Sj 
respectively denote the energy and entropy formation 
of the vacancy in pure A. 

The impurity diffusion coefficient, Dg, depends on 
several jump frequencies, corresponding to the ex¬ 
changes of the vacancy with the solute atom S and 
with the solvent atoms A near S, see Fig. [TJ 

D*s = o^^2fsCv , (3) 

where W 2 is the ^'-vacancy exchange frequency and fs 
is the solute correlation factor. 

For bcc lattices, in the formalism of Le ClaireQ, 
the correlation factor is 


jbcc 


1 -1 

TTt’ 


(4) 


where t is expressed in terms of the jump frequencies 
as: 


t= - 

i 02 T 3 cC 3 -f ScCg -f Ouj'f — 

with F = 0.512 i n (151) . 
correlation factor is 


—0J2 
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(5) 

For fee lattices, the solute 


ft = 
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( 6 ) 


with expressed as a function of it = a; 4 /a;o as 
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whith the coefficients calculated by Koiwa in [T^ . 

According to the transition-state theory, in a sys¬ 
tem of N atoms, the exchange frequency between a 
vacancy and an atom is. 


LOi = Vo exp 


ksTj 


= Vo exp 


TSm - K 
keT 


( 8 ) 



FIG. 1: (Color online) The frequencies involved in the sec¬ 
ond binding model for bcc and fee lattices. In black/orange 
circles, respectively are represented the solvent and solute 
atoms. 


In (IH) , G^i is the migration Gibbs free energy and the 
pre-exponential term, the "attempt frequency" vo, is 
of the order of the Debye frequency. The Gibbs free 
energy is given by Gm = Hm — TSm, where Sm is the 
migration entropy, while Elm is the enthalpy. As the 
volume is kept constant and the pressure is considered 
null, Hm = Em, where Em is the internal migration 
energy. Hence, following Vineyard’s formulation [l^ . 
the migration frequency jumps are given by 

w, = Vo exp(-£;;„/fcsT). (9) 

In ([S]) , El^ are the vacancy migration energies at T = 
OK, while 



with vl and vf the frequencies of the normal vibra¬ 
tional modes at the initial and saddle points, respec¬ 
tively. 

We present our numerical results applied to Fe- 
Cr diluted alloys. Above the melting temperature 
= 1183Ar, Fe-Cr alloys develop a paramagnetic 
fee phase, while for lower temperature the structure 
is bcc. In this bcc phase, a magnetic transitions oc¬ 
curs from ferromagnetic, below the Curie temperature 
Tc = 1043Ar, to paramagnetic states. 

In the case of the bcc phase, we performed both 
DFT and CMS calculations. For DFT calculations, 
we have employed localized basis sets as implemented 
in SIESTA code [l^. We have also considered spin 
polarization and GGA approximation in all calcula¬ 
tions. Core electrons are replaced by nonlocal norm- 
conserving pseudo potentials as in Ref. [l^. Valence 
electrons are described by linear combinations of lo¬ 
calized pseudoatomic orbitals. The basis sets for both 
elements consist in two and three localized functions 
for the 4s and 4p states, respectively, and five for the 
the 3d states. The maximum cutoff radius is 5.1 A. 
Calculations were carried out with 54 and 128 atom 
supercells, using respectively a 7 x 7 x 7 and 4x4x4 
k-point grid, and the Methfessel-Paxton broadening 
scheme with a 0.3eV width. The migration barriers 
have been determined using SIESTA coupled to the 
Monomer 0. 

In GMS calculations the atomic interaction are rep¬ 
resented by EAM potentials. For the Fe-Cr system in 






































3 


the bcc lattice we have used the potential developed 
by Mendelev et al. [2^, for the pure elements Fe and 
Cr, as well as, for the cross Fe-Cr term. While for the 
hight temperature fee phase, where only CMS calcu¬ 
lations were performed, we have used the potential 
developed by Bonny et al. [2l| . For all classical calcu¬ 
lations we use a christallyte of 8 x 8 x 8 with periodic 
boundary conditions, that is 1024 and 2048 atoms for 
bcc and fee respectively. We have verified, for the bcc 
phase, that the results do not change if we employ 
a christallyte of 128 atoms. We obtain the equilib¬ 
rium positions of the atoms by relaxing the structure 
via the conjugate gradients technique. The lattice pa¬ 
rameters that minimize the crystal structure energy 
are ape = 2.866A, and ape = 3.562A, for bcc and fee 
structures, respectively. 

In Table U we show our calculations of the activation 
energies, formation and migration, in a perfect bcc Fe 
lattice. We show both, DTF (with 54 and 128 atoms) 
together with CMS calculations. Initial and saddle 
points configurations and their respective energies are 
calculated with the Monomer method 0- We can 

TABLE I; Energies and lattice parameters for the pure bcc 
Fe lattice obtained by DFT calculations with 54 and 128 
atoms and by CMS calculations. 


bcc -Fe 


DFT 54 

DET\ 2 s 

CMS 

Exp. 

a{A) 

2.885 

2.885 

2.866 

2.866 

EjieV) 

2.18 

2.05 

1.72 

1.79 ±0.1 

E?nieV) 

0.67 

0.68 

0.68 


ElieV) 

0.57 

0.56 

0.562 


EUeV) 

0.67 

0.67 

0.67 


EUeV) 

0.64 

0.63 

0.625 


Et{eV) 

0.63 

0.60 

0.558 


Ei{eV) 

0.61 

0.60 

0.599 


E^ieV) 

0.60 

0.58 

0.542 


EtnieV) 

0.59 

0.59 

0.585 


ELieV) 

0.64 

0.63 

0.627 



observe from Table U that, as the number of atoms is 
increased, the energies calculated with DFT get closer 
to the classical ones. This effect is particularly impor¬ 
tant for the vacancy formation energy where the result 
obtained with CMS calculation is in accordance with 
the experimental result measured in (^ . 

For the fee paramagnetic phase, occuring above 
Ta^ = 11831^, formation and activation energies from 
CMS calculations in are displayed in Table HIl 


TABLE IL Activation energies in paramagnetic fee Fe-Cr 
from CMS calculations, using potential of Ref.. 


E^ 

Et 

Et 

Et 

Et 

Et 

Et 

■pA' 

Et 

^ 4 " 

1.87 

0.64 

0.66 

0.65 

0.76 

0.72 

0,66 

0.60 

0.70 

0.63 


In order to compute the jump frequencies, we use 
expression ®, with the migration energies E"!^ re¬ 
ported in Tables U and im respectively for bcc and 
fee phases. For the pre-factor we have taken the 

experimental values of the migration entropy Sm = 
2.1fcB, as reported in in all cases. While the De¬ 
bye frequency has been taken as vp = lO^^iLz. 


For the vacancy concentration in ([5]) , the formation 
entropy has been taken as S'J = d.lfc^, from DFT 
calculations performed in (l^l, while for CMS calcula¬ 
tions, we used Sj = 2.3fcB, as obtained in (2^. 

Once the jump frequencies in the multi-frequency 
model have been computed, the diffusion coefficients 
are calculated using analytical expressions © and ([3]). 
Also, it has been observed that, due to spontaneous 
magnetization [2^, the self-diffusion coefficient devi¬ 
ates from a linear Arrhenius relationship, below the 
Curie temperature. This magnetization effects are, as 
usually, taken into account, as a correction of the ac¬ 
tivation energies Q for the ferromagnetic phase, from 
those in the paramagnetic Qp (in Table |T]) such that, 

Q = Qp{lFaxs\{T)), (11) 

with X = Fe or Cr and sx {T) is the ratio of the s pon - 
taneous magnetization at T to that at T = OAT |^ . 
While Sx{T) = 0, in the full temperature range of the 
paramagnetic phase. In this respect, a direct estima¬ 
tion of these parameters from fir st p rinciples would 
be of great interest. Here, as in [^, we interpolate 
the values of as^{T) in Ref. for both, solute and 
solvent atoms. 

In Figured! we show the calculated D*p^ and D'^^, 
using equations o and respectively, with the ac¬ 
tivations energies in Tables HI and |ll| for bcc and fee 
phases, respectively. As we already mentioned, for 
the bcc phase we performed DFT, with 54 and 128 
atoms, as well as, CMS calculations. Also in figure 
dl experimental data taken from Refs. (2^ and [1^ 
are plotted with triangles and stars, respectively for 
D*p^ and As can be observer in Fig. dl below 



FIG. 2: (Color online) Self-diffusion (in black) and Cr im¬ 
purity (in orange) diffusion coefficients in Fe from CMS 
and DFT calculations. Full lines correspond to CMS cal¬ 
culations while dotted and dashed lines for DFT with 128 
and 54 atoms respectively. Experimental values for Dp^ 
and obtained from Refs. and are plotted 

with triangles and stars respectively. 

the solvent melting temperature T^-y = 1043°C', in ac¬ 
cordance with Bohr’s correspondence principle, as the 
size of the atomic cell (total number of atoms) is in¬ 
creased, quantum results with DFT approach to the 
classical ones obtained with CMS. 
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Note that the bee supercell with 54 atoms, corre¬ 
sponding to a christallyte of 3 x 3 x 3, has a length 
of 8 . 59 A, and for 128 atoms, the 4x4x4 christal¬ 
lyte has a length of 10.14A. In both cases, this length 
is lower than typical electronic quantum coherence 
leng th which are of nanometer order (4nm on Cm 
( 29|). Then, quantum coherence effects, which are 
only taken into account by DFT, play a crucial role 
for both the convergence issue of the DFT results as a 
function of supercell size, as well as, in the difference 
arising between classical CMS and quantum DFT cal¬ 
culations. 

Instead, with CMS, the semimpirical potential are 
phenomenological and the quantum coherence effects 
are not taken into account. Moreover, CMS calcula¬ 
tions with 128 or 1024 atoms, give the same results. 

In addition, we must not expect to observe quan¬ 
tum effects for such macroscopic systems especially at 
the high temperatures here described. In that case, 
the interaction of the system with the thermal en¬ 
vironment implies in decoherence effects, where the 
classical limit is expected to be recovered [1|. 

Indeed, our results obtained with CMS calculations 
are in good agreement with available experimental 
data for both, tracer solute and solvent diffusion coef¬ 
ficients. Note that, similar results have been obtained 
using, the also classical method. Kinetic Monte Carlo 
algorithm with temperature dependent pair interac¬ 
tions (30| . 

It must be emphasized that the agreement between 
CMS based calculations and experimentally measured 
diffusion coefficients is not fortuitous, it has been re¬ 
cently observed for diffusion in Al-U and Ni-Al fee 
lattices [l^- While for this former DFT calculation 
underestimated the diffusion coefficients 0. 

Several possible explanations of the fact that DFT 
calculation are not in agreement with ex^riments for 
the diffusion coefficients where argued in . We claim 


here that this is due to quantum coherence effects aris¬ 
ing from DFT calculations for the size of the simula¬ 
tion cell being small. As the size of the simulation cell 
is increased the DFT results converge to the experi¬ 
mental values that can be obtained with CMS calcu¬ 
lations which is much less expensive. 

For the fee phase, where diffusion coefficients have 
not yet been measured, our CMS calculations predicts 
the diffusion behavior. 

In summary, in this work we have performed a com¬ 
parison between quantum DFT and CMS calculations 
in order to obtain the diffusion properties in bcc Fe-Cr 
diluted alloys. In accordance with Bohr’s correspon¬ 
dence principle, as the total number of atoms is in¬ 
creased, the diffusion coefficients obtained with quan¬ 
tum DFT calculations, approach the classical ones ob¬ 
tained with CMS. For DFT calculations, the electronic 
quantum coherence plays a crucial role that is related 
with the size of the simulation cell. Also, thermal con¬ 
tact with the environment as the effect of killing coher¬ 
ence effects making the classical behavior to emerge. 
Indeed, results obtained with CMS calculations are in 
good agreement with available experimental data for 
both solute and solvent diffusion coefficients. 

Hence, the atomic diffusion process in metals is a 
classical phenomena, for which the large number of 
atoms and the temperature has suppressed any quan¬ 
tum coherent effect. Then, if reliable semi empirical 
potentials are available, a classical treatment of the 
atomic transport in metals is much convenient than 
DFT. 

The comparison between DFT and CMS calculation 
is then purposed as a tool to investigate the effective 
size of quantum effects, such as coherence length. 
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